Doppler Tomography of Relativistic Accretion Disks 



Vladimir KarasQ and Pavel KRAUS 



Astronomical Institute, Charles University, Svedskd 8, 
CZ-150 00 Prague, Czech Republic 



Scheduled for 
Publications of the Astronomical Society of Japan, 
Vol. 48 (October 1996) 



*Also at Scuola Internazionale Superiore di Studi Avanzati, Trieste; Department of 
Astronomy and Astrophysics, Goteborg University and Chalmers University of Technology, 
Goteborg 



Abstract 



Spectral lines from a source orbiting around a compact object are stud- 
ied. Time variations of observed frequency and count rate due to motion of 
the source and gravitational lensing are considered. Gravitational field of the 
central object is described by the Kerr metric. It is shown that: (i) simul- 
taneous temporal and frequency resolution enables us to restrict parameters 
of the model (inclination angle, position of the source, angular momentum 
of the black hole); (ii) techniques of image restoration, familiar from other 
fields of astronomy, can be applied to study inner regions of active galactic 
nuclei. This contribution is relevant for extremely variable X-ray sources with 
high parameter of efficiency, such as Seyfert 1 galaxy PHL 1092 observed by 
ROSAT. 
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1. Introduction 

Tomography is a method of image restoration which has been success- 
fully applied in various fields, particularly in medicine. Image reconstruction, 
on the other hand, is one of the typical problems in astronomy where the 
situation is rather specific because astronomical objects are very distant and 
can be viewed from a single direction only. Rotation of the source is thus 
essential, and one needs to understand how rotational motion is linked to 
observed profiles of spectral lines. Vogt, Penrod & Hatzes (1987) describe 
the tomographical method for mapping the surface of rotating stars. It is 
assumed that irregularities (spots) on the stellar surface radiate in a spec- 
tral line, observed frequency of which is periodically affected by the Doppler 
effect. An analogous approach was applied in order to study accretion disks 
in cataclysmic variable stars (Marsh, Home 1988; Kaitchuck et al. 1994) and 
Algol-type binaries (Albright, Richards 1996). It has been demonstrated that 
Doppler tomography can be very useful for tracing the structure of the disk 
surface when it is seen at a sufficiently large inclination. Similar techniques 
have been used recently to study detached binaries and multiple stars (Si- 
mon, Sturm 1994; Hadrava 1995). All these methods are based on the fact 
that the line spectrum of the source is affected by the radial velocity of its 
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individual components with respect to the observer. Knowing the compos- 
ite spectrum at different phases, the problem of the line formation can be 
inverted and contributions from separate components of the source can be 
revealed. 

The present paper deals with two topics: (i) the applicability of Doppler 
tomography to study rotation and surface structure of disks around compact 
objects is discussed, and (ii) the relevance of time-resolved spectral-line pro- 
files for study of the rapid X-ray variability of active galactic nuclei (AGN) 
is explored. (As for Doppler tomography, we retain the original terminology, 
although it is less appropriate in the context of general relativity in which the 
distinction between Doppler shift and the gravitational shift of the photon 
energy cannot be defined unambiguously.) The motivation for investigating 
these problems comes from a widely accepted model of AGNs with a central 
black hole surrounded by an accretion disk (see, e.g., Rees 1984). It has been 
realized by several authors that double-peaked asymmetric profiles of the Ha 
line observed in some galactic nuclei can be explained by relativistic effects 
(Chen, Halpern, Filippenko 1989; Eracleous, Halpern 1994). Theoretical line 
profiles have been studied by numerous authors who also discussed specific 
effects of general relativity: frame dragging (Kojima 1991; Laor 1991; Brom- 
ley, Chen, Miller 1996), multiple images (Bao, Hadrava, Gstgaard 1994), and 
self-gravity of relativistic disks (Karas, Lanza, Vokrouhlicky 1995). 

For a long time, some AGNs have been recognized as rapidly variable 
X-ray sources (Feigelson et al. 1986; Mushotzky, Done, Pounds 1993; Pa- 
padakis, Lawrence 1995). Apart from blazars, the class of I Zw 1 objects 
(Phillips 1976; Boiler et al. 1993) has been distinguished by the high value 
of their efficiency parameter: r] m 5 x 1CT 43 AL X /At <: 0.1 (Fabian 1979). 
This parameter characterizes temporal change of the source X-ray luminos- 
ity, ALx [ergs -1 ], during the relevant variability time-scale At « 10 2 -10 5 
seconds. For example, a member of the I Zw 1 class, narrow-line Seyfert 1 
galaxy PHL 1092 has recently been studied in detail (Forster, Halpern 1996) 
and it has been speculated that the remarkably high value of its efficiency 
parameter (rj > 0.6) calls for a non-standard model of the object (strongly 
anisotropic emission or accretion onto an extremely rotating black hole; Kwan 
et al. 1995). It is suggested in the present contribution that the lensing effect 
acting on a source of radiation orbiting a black hole (as discussed by numer- 
ous authors in another context; cf. Abramowicz, Bao, Lanza, Zhang 1991; 
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Rauch, Blandford 1994) should also be considered as a possible explanation 
of such extreme variability. The model can be submitted to an observational 
test by studying time-resolved spectral lines which are formed in that region 
of the source responsible for variability. Observations have not yet provided 
sufficient resolution for serious spectroscopy of this type. 

The most prominent relativistic effects can be expected for spectral 
features around 7 keV because this complex is presumably formed in the 
innermost regions of accretion disks (Fabian, George 1991; Matt, Perola, 
Piro, Stella 1992). The improved technology of ASCA's detectors makes high- 
resolution spectroscopy of the X-ray sources feasible (Serlemitsos et al. 1995) 
so that the above-mentioned theoretical works on line profiles now become 
timely. Possible future X-ray spectroscopy missions will provide temporal 
resolution of the spectral lines which has not yet been much discussed in the 
literature. The subsequent section of this paper deals with predicted line 
profiles as a function of time. Specific features of the frequency shift relevant 
for the tomographical method of the disk mapping are mentioned. The main 
aim of this discussion is to point out the differences between accretion disks 
around compact objects and other situations where the effects of relativity 
are negligible. 



2. Doppler Tomography of Accretion Disks 



2. 1 Details of the Model 

It is assumed that a localized region of enhanced emissivity is located on 
the surface of an accretion disk. The possible existence of such irregularities 
has been discussed by several authors because it has important consequences 
for variability of AGNs (Abramowicz, Lanza, Spiegel, Szuszkiewicz 1992; 
Sivron, Caditz, Tsuruta 1996). It is assumed that the disk resides in the 
equatorial plane of a Kerr black hole. It will be shown that the observed flux 
and frequency of a spectral line as functions of time determine, to a certain 
degree, the inclination of the disk and the angular-momentum parameter of 
the black hole, a (0 < a < 1). The relative frequency shift g of radiation 
(ratio of frequency observed by a distant observer to frequency emitted in a 
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local frame corotating with the disk material) can be expressed, within the 
approximation of geometrical optics, in terms of four-momentum of photons, 
p, and four- velocities u em , u ohs of the emitting material and of the distant 
observer, respectively: 

^obs P ' U obs / 1 \ 

9 = — = ; (i) 

^em P ' "em 

g is a function of polar coordinates {r, 0} in the disk plane, defined in such a 
way that {(f) = 0,9 = 9 ohs } represents a direction to the observer (9 ohs = 90° 
when the disk is seen edge-on). In order to minimize the number of free 
parameters, Keplerian rotation law is assumed: Q oc 1/ (r 3 / 2 + a) in usual 
dimensionless geometrized units. This restriction is in agreement with our 
assumption of geometrically thin, equatorial disk. The method can be applied 
to a different rotation law analogously. 

Radiation flux from the source is obtained by integrating the observed 
intensity J b s over the observer's local sky. The resulting value of the flux is 
determined by distribution of g and 7 e m (locally emitted intensity) in the disk 
(-^obs/-^em = 9 3 ) , and it also depends on the shape of light trajectories which 
are affected by the presence of the central object (lensing effect). Various 
semi-analytical and completely numerical approaches have been developed 
for efficient integration of the observed radiation flux (Pineault, Roeder 1977; 
Luminet 1979; Asaoka 1989). We employed the method described by Karas, 
Vokrouhlicky & Polnarev (1992). 



2.2 Line Profiles Resolved in Time 

First, it is instructive to approximate light rays by straight lines because 
the problem of observed line profiles can then be solved analytically without 
difficulty (Gerbal, Pelat 1981; Kojima 1991). Let us consider a point-like 
source moving along a circular orbit with radius r = r cm and emitting, in its 
rest frame, a spectral line with Dirac function-type profile: 

/ cm (z/,r,0) oc 5(v - v cm )5(r - r cm )5((f) - ilt). (2) 

The observed flux from a single source is a simple periodic function of time, 

F(u,t) oc g 3 I cm (v/g,r,(f) - r cos (f) sin 6 ohs ), (3) 
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and its observed phase ip can thus be normalized to the unit interval, < <p < 
1. The constant of proportionality in eq. (|]) is undetermined unless a specific 
physical mechanism generating radiation is assumed, but for the ratio of the 
count rates at the maximum and the minimum observed frequency one can 
write 

-^X^max) _ ^max ^max _ r cm Sm ^obs _ 
F^mm) ^min ^min T^m — SHI o bs 

r em can be estimated from relation (f|). Superposition of radiation from 
multiple sources results in a fluctuating signal but, in principle, individual 
contributions can still be restored. In practice, however, temporal and fre- 
quency resolution of observational data is limited, and one has to work with 
a more complicated local profile of the line. An appropriate technique of im- 
age restoration has been discussed by Marsh & Home (1988). These authors 
ignored any relativistic corrections (frequency shift was determined solely by 
radial velocity of the source), nevertheless, their approach can be applied to 
our situation if eqs. (|])-(@) for the redshift factor and intensity are taken into 
account. The difference becomes evident upon comparison of our Fig. [Dwith 
Fig. 1 of Marsh & Home: Notice the asymmetry of isocontours of g with 
respect to the vertical axis. The shape of the contour-lines is determined 
by mutual competition of the overall gravitational redshift with the Doppler 
effect due to rotation (right-hand side of the disk approaches the observer 
and the corresponding observed frequency is thus increased). As a result of 
this asymmetry, observed profiles of spectral lines are also asymmetric about 
the emitted frequency. 

The lensing effect, neglected in the discussion above, calls for a correc- 
tion of the observed flux. The magnitude of the error introduced by this 
neglect was discussed by several authors (Kojima [1991] corrects an incon- 
sistency in the work of Gerbal & Pelat [1981], but, in most situations, the 
difference does not produce a visible change in line profiles). Apparently, 
the correction becomes important when most of the radiation is generated 
near a black-hole horizon (typically, r < 20 r g , r g = 1 + yl — a 2 ) or when 
observer inclination is large (6* b s ^ 75°) . The approximation of straight pho- 
ton trajectories is unable to reflect lensing effects and we thus turn to a more 
sophisticated treatment by solving the geodesic equation in the Kerr metric. 

Fig. |2| illustrates variations of the line profile of a source in the Keplerian 
circular motion around a Schwarzschild black hole (a = 0). It is assumed 
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that the source has the form of a spot radiating isotropically in its local 
frame. Local emissivity has a Gaussian frequency profile (FWHM 0.2). 
Local intensity decreases exponentially with the distance from the centre of 
the spot. The characteristic diameter d of the spot is defined by requiring 
that intensity at the edge of the spot decreases by a factor 1/e in comparison 
with the value in the centre. In our example, d ~ 1 r g . The resulting observed 
profile is not sensitive to the assumptions above, provided 

2vrr cm »rf, FWHM« (5) 

^max ^min 

The local profile of the line is centred around z/ em but the observed profile 
oscillates around a lower value of frequency (observed centroid frequency) 
because of the gravitational redshift. The maximum frequency corresponds 
to the source approaching the observer (jp> « |7r), the minimum frequency 
corresponds to the source receding from the observer (<p rs ~7r). The observed 
count rate follows this variation strictly only if lensing and time delay are 
ignored (cf. eq. [[|). On the other hand, once these effect are take into 
account properly, the count rate has its maximum somewhere within the 
interval n < ip < |7r, depending on obs . Two well-separated local maxima 
of the count rate can be found in the signal. Lensing effect is responsible 
for enhancement at <p ~ 0.5 but the corresponding peak is visible only if 
#obs ~ 75° (Fig. ||a). Doppler shift produces a peak at <p « 0.7. One can 
estimate inclination by fitting the line profiles resolved in frequency and 
time, while with frequency integrated light curves such an estimate is rather 
uncertain (it can be shown that the phase of the maximum of frequency 
integrated signal depends only weakly on a and o b s ; cf. Fig. 8 of Karas, 
Vokrouhlick & Polnarev 1992). Since the geodesic equation was solved for 
each of the rays, the correct relation between the orbital phase <ft of the source 
on the disk surface and the the phase (p of the observed signal is automatically 
taken into account. In the case of a rotating black hole, line profiles take 
a similar course but the frequency range and phase of the maximum are 
different, as illustrated in the next figure. 

A relation between the observed frequency and the count rate is further 
explored in Fig. ||| where the centroid frequency of the line is given in terms 
of the count rate. The variable part of the count rate is normalized to unit 
interval. The centroid frequency is normalized to the local centroid frequency 
of the line profile in the frame comoving with the source. It is worth noticing 
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that the maximum and the minimum centroid energy (z/ Cjmax , f c ,mm) and the 
observed frequency of the maximum count rate (v) are specific indicators 
which relate each of the curves to the corresponding values of 9 em , r em , and a 
(Table 1). One can define the relative difference of the centroid frequencies: 

Av c = - 1 (6) 

^c,min 

which increases with # Q b s increasing and r cm decreasing. Large values of 
Ai/ C (> 0.5) correspond to relativistic motion of the orbiting source, while 
small V (< f Cjma x) indicate the focussing effect to be important. One can 
verify from Tab. f that, for large inclination, V decreases with r cm increasing, 
while for small inclination the dependence is reversed. Parameters Ai/ C and 
V supplement the asymmetry parameter A which was introduced by Gerbal 
& Pelat (1981, eq. [13]) in order to characterize asymmetry of the time- 
integrated line profiles. 



3. Conclusion 

The sufficient resolution, both in time and frequency, is crucial for prac- 
tical applicability of Doppler tomography and this is presently the main dif- 
ficulty of the method in studying inner regions of accretion disks in AGNs. 
In the present contribution we concentrated on properties of radiation from 
a single orbiting source, taking into account gravitational lensing of a Kerr 
black hole. This is the point where the situation differs from previous works- 
otherwise the techniques of image restoration can be applied in a similar 
manner as has been done for tracing, e.g., stellar surfaces or detached binary 
systems. 

Orbital motion of the source around a black hole results in variations of 
the observed count rate and increases the efficiency parameter of the source. 
This may be a relevant explanation for some extremely variable X-ray objects 
from the I Zw 1 class, and its relevance can be tested by studying line profiles 
in time. The method restricts possible values of inclination angle of the disk 
and angular momentum of the black hole, and it should thus be taken into 
account in planning future X-ray observations. 
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We assumed that the source of light radiates isotropically in its local 
comoving frame, and that it rotates in accordance with the Kepler law. The 
assumption about local emissivity does not play a particular role in our 
discussion of the line profiles (Fig. |l|) because we did not consider general 
relativistic bending of light trajectories. Fully relativistic treatment of Figs. 

is however affected by this constrain. Also the rotation law needs to 
be determined in accordance with the physical situation (equatorial motion 
around a Kerr black hole, in our case). Otherwise, it introduces another 
free parameter which must be determined by systematic exploration of the 
parameter space and fitting values from Tab. 1 to actual data. 

This work has been supported by the grants GACR 205/94/0504 and 
GACR 202/96/0206 in Czech Republic, and by Wenner-Gren Center Foun- 
dation in Sweden. 
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Figure Captions 



Fig. m Lower panels: Lines of equal frequency shift g on an equatorial accre- 
tion disk in Keplerian rotation around a Schwarzschild black hole. Observer 
looks from the bottom of the page (0 = 0). Both graphs are drawn with 
radial coordinate x = 1 — 3r g /r (r g = 2, < x < 1). The whole equatorial 
plane outside the inner edge of the disk is thus captured in this illustration; 
dotted circles corresponds to x = 1. Observer inclination is (a) 80° and (b) 
20°. Upper panels: Observed count rates (eq. 0) from a point-like source 
orbiting at r em = 3r g (solid line) and r em = 8r g (dashed line). The count 
rate is in arbitrary units; frequency is normalized to z/ em . 

Fig. |2|. Time variation of the line profile from a source near a Schwarzschild 
black hole; (a) r em = Ur g , 8 ohs = 80°; (b) r em = 14 r g , 6 ohs = 20°; (c) 
r em = 4r 9 , (9 b s = 20°. Count rate is normalized to the maximum value; 
frequency is normalized to z/ em . The phase interval corresponding to one 
revolution of the source is shown (0 < ip < 1). See the text for details. 

Fig. |3|. Centroid frequency versus count rate. Each loop characterizes one 
complete revolution of the source at given r em . Three different radii of the 
orbit have been chosen, as indicated by the line-style: r em = 3r g (solid 
line); r em = 5r g (dashed line); r em = 8r g (dot-dashed line). Four panels are 
shown: (a) 9 ohs = 80°, a = 0; (b) # obs = 20°, a = 0; (c) # obs = 80°, a = 1; (d) 
#obs = 20°, a = 1. Results from this figure are summarized in Tab. 1. 



Hardcopy of figures are available upon request from the authors. You can 
also retrieve the complete postscript file of this paper. 



Figure 1: V. Karas & P. Kraus 



Figure 2: V. Karas & P. Kraus 
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Figure 3: V. Karas & P. Kraus 
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Table 1. Characteristics of the centroid frequency 
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a = 









a = 


1 




^c,max 


^c,min 


Au c 


V 


^c,max 


^c,min 


Au c 


V 


80 


3 


1.32 


0.53 


1.52 


1.13 


1.31 


0.51 


1.55 


1.19 


20 


3 


0.84 


0.60 


0.39 


0.83 


0.68 


0.52 


0.29 


0.67 


80 


5 


1.26 


0.62 


1.05 


1.00 


1.31 


0.51 


1.55 


1.08 


20 


5 


0.95 


0.74 


0.28 


0.95 


0.83 


0.59 


0.41 


0.83 


80 


8 


1.23 


0.71 


0.74 


0.99 


1.27 


0.58 


1.20 


0.93 


20 


8 


0.98 


0.82 


0.18 


0.98 


0.92 


0.71 


0.28 


0.92 



